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Abstract. We investigate resonant strong field ionization of atomic hydrogen with respect 
to the ls-2p-transition. By "strong" we understand that Rabi-periods are executed on a 
femtosecond time scale. Ionization and AC Stark shifts modify the bound state dynamics 
severely, leading to nonperturbative signatures in the photoelectron spectra. We introduce an 
analytical model, capable of predicting qualitative features in the photoelectron spectra such 
as the positions of the Autler-Townes peaks for modest field strengths. Ab initio solutions of 
the time-dependent Schrodinger equation show a pronounced shift and broadening of the left 
Autler-Townes peak as the field strength is increased. The right peak remains rather narrow 
and shifts less. This result is analyzed and explained with the help of exact AC Stark shifts 
and ionization rates obtained from Floquet theory. Finally, it is demonstrated that in the case 
of finite pulses as short as 20 fs the Autler-Townes duplet can still be resolved. The fourth 
generation light sources under construction worldwide will provide bright, coherent radiation 
with photon energies ranging from a tenth of a meV up to tens of keV, hence covering the 
regime studied in the paper so that measurements of nonperturbative, relative AC Stark shifts 
should become feasible with these new light sources. 



1. Introduction 

Resonant two-photon ionization of atomic hydrogen starting from the ground state has been 
studied for more than thirty years as a prime example for theoretical models in which two 
discrete states are coupled to the continuum (see, e.g. [[Ill2l[3l|4l|51|6l|7l[H an ^ references 
therein). The photoelectron spectra display duplets due to the AC Stark splitting of the bound- 
bound resonance (Autler-Townes [9J duplets). When plotted as a function of the detuning the 
energies of the two Autler-Townes peaks directly map the two field-dressed state energies into 
the continuum, thus making measurements of AC Stark shifts and avoided crossings possible 
ifTOl [TTfl . Besides this, it has been shown recently that the interference in Autler-Townes 
duplets can be controlled using two time-delayed intense femtosecond pulses lfl2l . 

If the coupling to the continuum is neglected, AC Stark shifts are ignored (or considered 
as input governing the effective detuning) and the rotating wave approximation is adopted, 
the well-known, analytically soluble Rabi-flopping dynamics are obtained. The influence of 
the latter on the emission spectra has been thoroughly studied for modest laser intensities (see 
Refs. |[T3l[T4| for reviews). More recently, harmonic generation from ionization-damped two- 
level systems [fl31 . radiation emitted by a resonantly driven H atom tfT6l and the enhancement 



Nonperturbative resonant strong field ionization of atomic hydrogen 



2 



of intense-field high harmonic generation due to the coherent superposition of states IfTTl [T8l 
have been studied theoretically. 

The unperturbed Rabi-flopping may be considered the zeroth order solution to be used 
for the calculation of spectra. Such approaches have been recently pursued to investigate 
ion impact ionization in the presence of a resonant laser [fT9ll and ionization of positronium 
in short, resonant UV laser pulses EOl . This perturbative approach will work as long as 
ionization and AC Stark effect do not influence the bound state dynamics significantly. Our 
current work aims at studying resonant strong field ionization beyond perturbation theory. On 
an ab initio basis and for short pulses this can be only achieved by a full numerical solution 
of the time-dependent Schrodinger equation (TDSE) ^2A\. 

Our work is inspired by the so-called fourth generation light sources now under 
construction worldwide (e.g. at DESY, Germany, in the UK and in the USA). These new light 
sources will provide coherent, short- wavelength radiation at intensities where the assumption 
of an unperturbed bound state flopping is invalid, opening up new possibilities to actually 
measure nonperturbative, relative AC Stark shifts of low-lying atomic or ionic states. First 
experiments on intense few-photon ionization of atoms have been already performed at the 
FLASH facility at DESY ||22ll23ll. 

The paper is organized as follows. In Sec.[2]we introduce an analytical model, capable of 
describing the main qualitative features in the photoelectron spectra, which are presented in 
Sec. [3] In Sec.|4]a& initio results from the solution of the TDSE are presented and compared 
with the model results. The nonperturbative features in the TDSE spectra are analyzed and 
explained with the help of Floquet theory in Sec. [51 The case of strong field resonant ionization 
by femtosecond pulses is studied in Sec.[6l Finally, we conclude in Sec. [71 



2. Model 



We consider a Hamiltonian of the form (atomic units are used unless specified otherwise) 

H(t) = H + W(t), H = T + V (1) 

where W(t) = zE(t) = zE(t) coscut accounts for the interaction with the linearly polarized, 
resonant (or almost resonant) laser field E(i) = E(t)e z in dipole approximation. We 
approximate the solution of the TDSE 

i|*(t)> =#(*)!*(*)> (2) 
by using the ansatz 

\m(t)) = e - i£at a{t)\a) + e- i£bt b(t)\b) + J d 3 k c{k,t)\k) . (3) 

Omission of the third term on the right hand side would lead us to the standard Rabi-theory 
(see any text book on quantum optics, e.g. ff24|). The ansatz the neglect of continuum- 
continuum transitions ~ (k'\V\k) (rescattering) and the assumption of mutual orthogonality 
of the states \a), \b), \k) is equivalent to the use of the model Hamiltonian 

H(t) =H + W(t) (4) 
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H = £ a \a)(a\+£ b \b)(b\+ Jd 3 k^-\k)(k\ (5) 

and 

W(t) = W hh (t) + W bc (t) + W cc (t) (6) 

where the subscripts "bb", "be" and "cc" indicate bound-bound, bound-continuum and 
continuum-continuum coupling due to the laser field E{t), respectively. The terms are 
explicitly given by 

W hh (t)=E(t)(a\b}(a\+h.c), 

W hc (t)=E(t) Jd 3 k(d a (k)\k)(a\+d b (k)\k)(b\ + h.c), 

W cc (t) = \E(t) Jd 3 kd 3 k' \k)d k J (3 \k - k')(k'\ 
with the transition matrix elements 

a = (b\z\a), d a (k) = (k\z\a), db(k) = (k\z\b). 

The model Hamiltonian TC(t) is hermitian so that for, e.g., the initial conditions a(0) = 1, 
6(0) = c(fc,0) = normalization |a(t)| 2 + \b(t)\ 2 + J d 3 k \c(k,t)\ 2 = 1 is ensured for all 
times. 

The three coupled equations governing the time-evolution of the coefficients a(t), b(t) 
and c(k, t) read 

id(t) = E(t) (a*b(t)e-^ Sb - £a)t + e i£at Jd 3 k c{k, t)d* a {k)^j (7) 

ib(t) = E(t) (aa{t)e i(£b - £ ^ + e i£bt Jd 3 k c{k, t)d* b {k)^j (8) 
k 2 

ic(fe, t) = —c{k, t) + K (a, 6, k, t) + iE(t)d kz c(k, t) (9) 



where 



-iS h t 



K(a, b, k, t) = E(t) (a(t)d a {k)e~ i£at + b{t)d b {k)e~ 

The partial differential equation © can be formally integrated using the method of 
characteristics. The result reads 

c(k, t) = -i f dt' K(a, b,k + A{t') - A{t), t')e- iS *(*'*'> (10) 

with A(t) = - f dt' E(t') the vector potential and S k (t, if) the action, 

S h (t,t') = \ fdt" [A(t") - A(t) + k] 2 . (11) 

Introducing X = (ae~ l£at , be~ 1 ( £b ~ u; ' t ) T , assuming a = a* being real and applying the 
rotating wave approximation (RWA) to the bound-bound dynamics lead to 



iX 



( p 9* \ 

VT~ £h ~ u J 
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where 

n R = Ea } (13) 
I a (t)= Jd 3 kd* a (k)c(k,t), (14) 

I b (t) = J d 3 kd* b (k)c(k,t). (15) 

If the second term on the right hand side of Eq. (fT2j) were known as an explicit function of 
time the differential equation could be solved analytically. However, since c(k, t) depends via 
n(a, b, k, t) on a(t) and b(t), Eq. (fT2j) is an integro differential equation still too complicated 
to be solved analytically. Close to resonance and for modest ionization one may neglect in 
zeroth order the second term, obtaining the well-known unperturbed Rabi-flopping dynamics 
(see, e.g. El) 



a (t) 



cos(fit/2) + ^sin(fit/2) 



e" iA */ 2 , (16) 



b (t) = ^ sin(fit/2)e iA '/ 2 (17) 
with 

n = yjnl + A 2 , A = £ b -£ a -u, (18) 

the Rabi-frequency and the detuning, respectively. We can then use (fT6l) . (fTTb for a(t) and b(t) 
in k(cl, b, k, t) to calculate c(k, t) in first order according (fTOl) by simple numerical integration 
over time (for each k of interest). Subsequently, the photoelectron spectra can be calculated 
from the momentum distribution P(k) = \c(k, t — > oo) | 2 . The vector potential can be chosen 
such that A(t — > oo) = 0. The probability for the emission of an electron with energy 
£ = k 2 /2 into the solid angle element df^ is then given by 

W - S§ - kPik ^ (19) 

In the weak-field limit one may neglect the vector potentials in the argument of k and in 
the action (fTTI) . Equation (flOl) then simplifies for t — >• oo and A = up to a constant phase to 
the first order perturbation theory (in uj) result 

c (k) = r&tfE{t') (a (t')d a (k)e- i£at ' + b (t')d b (k)e~ iSbt ') e ifcV/2 

~ E f°°dt' d b {ky^' 2 - £b -^ U^ 1 '/ 2 - e-^'l 2 ) (20) 

4 J — oo ^ ' 

where the subscripts '0' indicate that the unperturbed Rabi-dynamics governed by Eqs. (fT6b . 
(fTTT) have been used. From Eq. (l20l) follows that a pair of peaks (a so-called Autler-Townes 
duplet) in the photoelectron spectra is expected at 8 = 8 b + u ± ^r/2. Generalizing this 
result to higher photon orders we obtain 

S = S b + nu±^, n = l,2,.... (21) 
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3. Model results 



For the case of atomic hydrogen, interacting with a linearly polarized laser field that is 
resonant (or almost resonant) with the Is <-> 2p -transition, i.e. £& — £ a = 0.375 = u + A, 
we calculated Pf,{£) according Eqs. (fl~9l) using (flOl) with the populations (fl~6l) . (flTT) and the 
dipole matrix elements for hydrogen-like ions of charge Z (in our case Z = 1), 

d.(*) =*.(*) = -bW_£_, (22) 

The Rabi-frequency at resonance is given by 

256 E 



n 



R 



243 v^' 



(23) 



(24) 



32 

>- 
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Figure 1. Logarithmically scaled photoelectron spectra (fl9l l in polarization direction k = e z 
for i? = 0.04 at resonance (a; = 0.375). For simplicity, a rectangular 32-cycle pulse was 
assumed. The atom undergoes ~ 2.5 Rabi-cycles. The expected positions of the two Autler- 
Townes peaks (1211 are indicated by the dashed vertical lines, the center of them by the dotted 
vertical lines. 

Figure [H shows the spectrum at resonance (u = 0.375) for a rectangular 32-cycle pulse 
of E = 0.04. During such a pulse the atom undergoes ~ 2.5 Rabi-cycles. The positions of 
the Autler-Townes duplets, as expected from Eq. (T2T)) . are confirmed. 

It is instructive to trace the transition from the nonresonant to the resonant two-photon 
ionization. Figure [2] shows the photoelectron spectra for various detunings A between the 
limits of double and single photon ionization, i.e. uj = 0.25 and to = \£ a \ = 0.5 or 
A = ±0.125, respectively. The expected nonresonant lowest order peak positions follow 
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straight lines and are given by p a (u>) = £ a + 2u, Pbiyj) = £ b + uj. These lines cross 
at resonance while the two actual photoelectron peaks do not. This is an experimentally 
observable manifestation of the well-known avoided crossings of field-dressed states lfT0l[TTT| . 
The two peaks appear to "collide" with (and repel) each other. Their separation at the closest 
approach (i.e. at resonance) equals the Rabi frequency Or [cf. Eq. ([2T]) ]. 



0.10 

-J 0.05 
d 

■ - 0.00 

Q 

— 0.05 

— 0.10 



\ v\ 















































































































I I , , , , 


\ \ , 



0.0 0.1 0.2 0.3 0.4 0.5 

Photoelectron energy (a.u.) 



Figure 2. Spectra for various detunings A = 0.375 — u> (other parameters as in Fig. [TJ. 
The red, blue and green lines indicate £ a + 2uj, £\, + u> and £ a + 2w res , respectively, with 

W res = 0.375. 



The model introduced in Sec. [2] relies on several simplifying assumptions. By taking 
the two-level dynamics (fT6l) . (flTl) as the zeroth order solution to our problem of an infinite 
number of bound states plus a continuum coupled by a laser field and making the RWA, we 
neglect ionization losses, AC Stark shifts and anti-resonant terms. 

It is well-known that ionization (or other losses) may be phenomenologically 
incorporated in the dynamics of a two-level system by introducing complex energies and 
complex Rabi-frequencies (see, e.g. Il25l and references therein). In our model, this 
corresponds to the assumption that the second term on the right hand side of (fl"2l) can be 
rewritten in the form 



E(t) 



Ia(t) ' 



(, X a .fii \ 

' r — 1 1 



a, 



— db-i— 



X 



(25) 



2 / 



so that 



iX 



1 ' £ a + S a - i y 
V 2 



£ b — uj + 5b — i— 



X. 



(26) 



2 / 
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Here, 5 a ^ are the AC Stark shifts and T a b are the ionization rates. Because of 

f t [\a(t)\ 2 + \b(t)\ 2 ] = -T a \a(t)\ 2 -T b \b(t)\ 2 

+ 2nM[a*(t)b(t)e-' lAt ] 

the bound state population then overall decreases due to ionization, provided f2i is chosen 
sufficiently small. In Ref. E51 the parameters 5 a ,b, r a ,6 and f2; were determined from Floquet 
calculations for laser field strengths smaller than the ones of interest here. 

4. Results from the numerical ab initio solution of the TDSE 

The validity of the above model can be checked by comparing to ab initio solutions of the 
TDSE i|\l/) = H(t)\ty). We used the Qprop package E6ll to propagate the exact wave 
function and to calculate the photoelectron spectra. Figures [3] and |4] show the TDSE results 
corresponding to the results of Figs. Q] and [2l respectively. While the qualitative agreement 
is good as far as the peak positions are concerned, the relative strengths of the two peaks are 
different in the TDSE and the model calculations. Comparing Figs. [2] and |4l it is clearly 
seen that the strength of the peak p b is overestimated in the model for positive detuning 
uj < Eb — S a . This is in agreement with observations made in Ref. (|2T]| where also the 
left Autler-Townes peak was found to dominate the right peak. Moreover, the TDSE results 
show that ionization from the Is state prevails already for small positive detunings, justifying 
one of the assumptions of the so-called strong field approximation ll27ll where all bound states 
different from the initial state are neglected. 




Photoelectron energy (a.u.) 



Figure 3. TDSE electron spectra in polarization direction k = e z for E = 0.04 and 
uj = 0.375, to be compared with the model result in Fig. Q] The laser pulse was 32 cycles 
long and of trapezoidal shape (up and down ramping over one cycle). 
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Figure 4. TDSE spectra for various detunings A = 0.375 — u> (other parameters as in Fig. [3j, 
to be compared with the model result in Fig.|2] 
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Figure 5. Ground state populations |a(i)| 2 vs time (in cycles of the driving laser of resonant 
frequency w = 0.375) for E = 0.035, 0.076 and 0.17. 

Next, we investigate how the photoelectron spectra change when the atom is resonantly 
driven stronger and stronger, up to E = 0.17 (corresponding to an intensity ~ 10 15 Wcm~ 2 ). 
Here, by "resonant" we understand a constant laser frequency of o> = 0.375, which is the 
"true" resonance frequency only for negligible AC Stark effect and ionization, i.e. 5 a ^ ~ 
~ r a> 6. The AC Stark effect is expected to push the system towards a positive detuning 
A = £ b + 5b — {£ a + S a ) — oj > since 5b > 5 a . Hence, the dynamic AC Stark effect leads 
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Figure 6. Angle-integrated photoelectron spectra for various laser field amplitudes E (given 
to the right) for co = 0.375. The spectra were vertically shifted for better visibility. The 
Autler-Townes duplets i = 1,2,3 are indicated. 




Figure 7. (a) Peak positions of the first three Autler-Townes duplets vs electric field amplitude 
(u; = 0.375) as obtained from the TDSE solution (symbols) compared to the expected 
positions according £\, + nu ± £Ir/2 with n = 1, 2, 3 (dashed), (b) Relative peak separations 
(di — SlrO/f^R of the two Autler-Townes peaks in the zth Autler-Townes duplet for i = 1,2,3. 



to an increasing Rabi frequency. On the other hand, ionization losses damp the system, and 
damping is expected to shift the Rabi-flopping frequency to smaller values. 

Both AC Stark shifts of all states (inch the continuum) as well as ionization are included 
in the TDSE solution. In the following TDSE calculations a trapezoidal 100-cycle pulse 
was applied, ramped up and down over one cycle. For the lowest laser intensity shown 
(E = 0.035) about 6.5 Rabi-cycles occur during the pulse while at the highest intensity (E = 
0.17) ionization is already so violent that a Rabi-period cannot be identified unambiguously 
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Figure 8. Separation of the two peaks in the first Autler-Townes duplet (i.e. the "true" Rabi- 
frequency IT) as a function of the laser frequency for E = 0.04 (blue, +) and E = 0.12 
(red, *). The expected result according d27T i (dotted) and the shifted parabolas d29l ) (dashed) 
are included. From the latter, the "true" detuning A and resonant Rabi frequency Or can be 
determined. 

anymore. Figure [5] shows the ground state populations | a. (t) | 2 vs time for E = 0.035, 0.076 
and 0.17, respectively. While, as expected, the Rabi-frequency increases with increasing field 
strength, ionization leads to losses, i.e., the ground state population |a(t)| 2 does not return to 
unity. In the strongly over-damped case (bottom plot in Fig. [5]) ionization clearly dominates 
the bound state dynamics. As long as <C Or the oscillation frequency of the population 
transfer between the Is and the 2p state is in good agreement with the separation of the two 
peaks in an Autler-Townes duplet in the photoelectron spectra. We determined the separation 
dj of the two peaks within the i-th Autler-Townes duplet from the photoelectron spectra 
shown in Fig. [6l As expected, the separation increases with increasing laser field strength. 
Additional peaks originate from ionization involving higher- lying bound states. The left peak 
of each Autler-Townes duplet broadens as the field amplitude increases while the right Autler- 
Townes-peaks remain narrow. One may naively expect it to be the other way around because 
for increasing detuning the left (right) peak becomes the £ a + 2opeak (£ b + opeak) and 
ionization from the excited state \b) should be more probable, leading to a broadened right 
peak. Why the contrary is true will become clear in the next Section where we analyze the 
results in terms of Floquet theory. 

In Fig. U\ the (a) Autler-Townes peak positions and the (b) relative separations (di — 
f2 R )/^R of the two peaks in the i-th Autler-Townes duplet for i = 1,2,3 are plotted as 
a function of the driving laser field amplitude. With increasing driver strength the Autler- 
Townes duplet is shifted towards energies lower than those expected from the simple formula 
(|2T|) . Figure Ub shows that the peak separation di becomes smaller than ^r, indicating that 
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the frequency down-shift due to ionization dominates the frequency up-shift due to the AC 
Stark shift detuning. 

Figure [8] shows the separation d i= i of the two peaks in the first Autler-Townes duplet, i.e. 
the "true" Rabi-frequency ti, as a function of the laser frequency for E = 0.04 and E = 0.12. 
In the modest intensity case E = 0.04 (corresponding to 5.6 x 10 12 Wcm~ 2 ) the TDSE-result 
is close to the expected result according 

Q(uj) = + A 2 (uj) (27) 
(included dotted in Fig. [8]). A horizontal shift corresponding to the replacement 

A -> A(w) = S b + S b - {S a + 5 a ) - uj = AS - uj + 6 b - 6 a (28) 

with AS = 0.375 the unperturbed level spacing and 5 b — 5 a = 0.007 yields the good agreement 
with the dashed curve ti(u) = \Jq r + A 2 (uj) in Fig. [8] for E = 0.04. For a pronounced 
negative detuning uj > 0.375 + S b — 5 a resonances with other (higher- lying) excited states 
come into play. In fact, in the TDSE spectra Autler-Townes duplets corresponding to the 
resonant coupling of the Is ground state with the 3p (A£i s _ 3p = 0.4) emerge when uj is closer 
to 0.4 than to 0.375 + 5 b - 5 a . 

At higher laser intensities the Rabi frequency is shifted towards lower values owing to 
ionization (cf. Fig. 05). Hence, by fitting the TDSE results in Fig.[8]to 

ti(uj) = v^r + A 2 H (29) 

with A(lu) defined in (|28b we can determine both the relative AC Stark shift 5 b — S a and the 
"true" resonant Rabi-frequency fi R . In the strong field E = 0.12-case (5.1 x 10 14 Wcm -2 ) a 
clear minimum (i.e. resonant Rabi-frequency tin) and horizontal shift (i.e. A — A) cannot be 
determined anymore. 

5. Floquet results 

The TDSE results for long driving pulses may be best analyzed and understood in terms of 
the complex Floquet energies e = S + 5 — iT/2 which were introduced phenomenologically 
in Eq. (l26l) . We used the STRFLO code [28J to determine the exact Floquet energies. 

Figure [9] illustrates how 9?(e) = S + 5 and T = —23(e) (i.e. the ionization rates) behave 
as a function of the driver frequency for the modest field strength E = 0.04. Well below 
the expected resonance at uj = 0.375 the energies follow S a + 5 a and S b — uj + 5 b . The 
AC Stark shift of the Is state is negative while the 2p state is shifted upwards by ~ U p [the 
dashed-dotted line is S b — uj + 5 b with 5 b = U p (u)]. Resonance occurs at the frequency 
where the two energies lie closest together. Since 5 b — 5 a > the resonance is shifted from 
the unperturbed value uj = 0.375 to the higher value ~ 0.38. The line £ n =3 — ^ crosses 
S a = —0.5 at the ls^3p resonance at uj = 0.4. In Fig. [9}) the corresponding ionization rates 
do cross at uj ~ 0.35, i.e. already below resonance. For frequencies 0.25 < uj < 0.3 the 
rates are well approximated by log 10 T a = —5uj — 1.89 and log 10 T b = —Iuj — 1, respectively 
[straight lines in (b)]. For uj = 0.375 the ionization rate of the lower branch (drawn black 
Fig-Hk) exceeds the ionization rate of the higher-lying red branch. This explains why the left 
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peak in an Autler-Townes duplet is broader than the right peak although the left peak is the 
one belonging to ionization from the 1 s state as one detunes the frequency towards smaller 
values. However, close to resonance the dressed states are superpositions of the unperturbed 
states with sizable contributions from both \a) and \b). 
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Figure 9. Real parts of the relevant Floquet energies as a function of the driver frequency u> 
for E — 0.04 (a) and the corresponding ionization rates (b). 

Figure [TOl shows the Floquet energies and ionization rates for the stronger field E = 0.12. 
The AC Stark shifts are very pronounced, shifting the resonance further towards higher 
frequencies where, however, the ls^3p resonance clearly affects the Floquet energy of the 
upper branch (red) whose AC Stark shift 5 6 is less than U p for such high field strengths. 
With increasing frequency the branches approach Sb — to (red, dotted) and £ n=3 — u> (black, 
thin solid), respectively, with the mutual distance £ n=3 — £ b = 0.069 (indicated by the 
arrow) becoming constant. This explains why the peak separations obtained from the TDSE 
simulations shown in Fig. [8] both tend towards the value ~ 0.07 with increasing frequency. 
The ionization rate of the lower branch is higher than that of the upper branch in the whole 
frequency range shown, not just close to resonance. This is a manifestation of so-called 
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adiabatic stabilization ||29l which is known to set in for driving frequencies exceeding the 
ionization potential and for sufficiently high field strengths. The pronounced AC Stark down- 
shift of the lower branch and the high ionization rate associated with it are responsible for the 
respective shift towards lower energies and the broadening of the left Autler-Townes peak in 
the TDSE results of Fig. [6l An analysis of the Floquet wavefunctions corresponding to the 
energies shown in Fig.[lO]reveals that at such high field strengths more states than just \a) and 
| b) contribute in the frequency range of interest. 
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Figure 10. Same as in Fig.|9]but for E = 0.12. (a) With increasing frequency the branches 
approach £b — to (red, dotted) and £ n —3 — u> (black, thin solid), respectively, with the mutual 
distance £ n —3 — £b — 0.069 (indicated by the arrow) becoming constant, (b) The ionization 
rate of the lower branch is higher than that of the upper branch (adiabatic stabilization). 

We have checked that the use of the AC Stark-shifted energies and ionization rates 
obtained from Floquet theory in Eq. (l26l) reproduces well the exact evolution of the 
populations |a(i)| 2 and \b(t)\ 2 from the TDSE. However, the quantitative agreement between 
the TDSE spectra and the model spectra calculated from (flOl) using these exact populations 
is poor as far as the relative strengths of the Autler-Townes peaks and absolute figures for, 
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e.g., the ionization probability are concerned. Giving up the RWA does hardly change the 
model results and thus does not improve the agreement, meaning that Bloch-Siegert shifts are 
of minor importance. The reason for the quantitative disagreement between model and TDSE 
results is most probably the plane wave ansatz in Eq. © (instead of Coulomb continuum 
states). It is known that the same approximation plagues the strong field approximation, and 
attempts to include Coulomb effects in the latter have been pursued (see, e.g. Refs. 113011311 
and references therein). 
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Figure 11. Same as in Fig.|6]but for a finite pulse E(t) = E cxp[— 4(ln 2)i 2 /Tp WHM ] cosujt 
with T FWH m = 20 fs. 



6. Finite-pulse results 



The new fourth generation light sources will deliver bright, coherent, short wavelength 
radiation as short pulses. It is clear that in cases where no Rabi-floppings can be completed 
within the pulse duration no clearly separated Autler-Townes peaks can be expected. 
According Eq. (1241) Rabi frequencies on the time-scale of a few femtoseconds are expected 
already at laser intensities ~ 3.5 x 10 14 Wcm~ 2 so that in a, say, 20-fs pulse many Rabi- 
cycles are executed and a clear separation of the two Autler-Townes peaks is anticipated. 
However, ionization and the AC Stark effect lead to a pronounced broadening of the left 
Autler-Townes peak at high intensities even for a flat-top pulse, as is visible in Fig. |6] Hence 
the question arises whether in a real high-field, short-pulse experiment the Autler-Townes 
duplets can be resolved. Figure HH shows the TDSE results for exponential 20-fs FWHM 
pulses (with respect to the electric field). Twenty femtoseconds is a realistic pulse duration 
for the new short wavelength radiation sources. The same peak field strengths as in Fig. [6] 
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were used so that flat-top and finite-pulse TDSE results can be directly compared. As to 
be expected, a broadening and smoothing due to the now time-dependent Rabi-frequency 
JIr = Q-n[E(t)] is visible. Nevertheless, the two Autler-Townes peaks are clearly separable, 
the contrast becoming better for higher order Autler-Townes duplets, presumably because of 
Coulomb effects having less influence on the more energetic photoelectrons. 

7. Conclusions 

We studied strong field ionization of atomic hydrogen close to the ls-2p-resonance by means 
of an ab initio solution of the time-dependent Schrodinger equation. The photoelectron 
spectra directly reflect properties of the field dressed states, e.g. avoided crossings at 
resonances, different shifts of the two Autler-Townes peaks due to the highly nonperturbative 
AC Stark effect, broadening of peaks due to ionization and higher order Autler-Townes 
duplets. We found ionization rates and AC Stark shifts hardly accessible to analytical theory 
for field strengths where the Rabi-dynamics occur on a femtosecond time scale. However, 
we were able to explain the partly counter-intuitive results (such as the broadening of the left 
instead of the right Autler-Townes peak) in terms of complex Floquet energies, giving the 
exact AC Stark shifts and ionization rates for infinite laser pulses. Finally, we showed that the 
Autler-Townes duplets remain resolvable in the spectra even when femtosecond pulses are 
applied. This is important in view of the fourth generation light sources under construction 
worldwide. We propose to use these novel sources of bright, coherent radiation to investigate 
strong field resonant ionization of low-lying atomic or ionic states experimentally in the 
nonperturbative regime. 
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